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Abstract 



Single crystal films of bismuth-substituted ferrite garnets have been synthe- 
sized by the liquid phase epitaxy method where GGG substrates are dipped 
into the flux. The growth parameters are controlled to obtain films with 
in-plane magnetization and virtually no domain activity, which makes them 
excellently suited for magnetooptic imaging. The Faraday rotation spectra 
were measured across the visible range of wavelengths. To interprete the 
spectra we present a simple model based on the existence of two optical tran- 
sitions of diamagnetic character, one tetrahedral and one octahedral. We find 
excellent agreement between the model and our experimental results for pho- 
ton energies between 1.77 and 2.53 eV, corresponding to wavelengths between 
700 and 490 nm. It is shown that the Faraday rotation changes significantly 
with the amount of substituted gallium and bismuth. Furthermore, the ex- 
perimental results confirm that the magnetooptic response changes linearly 
with the bismuth substitution. 



I. INTRODUCTION 



Bismuth-substituted ferrite garnets (Bi:FGs) are well known to have a giant magnetoop- 
tical response For this reason, they have found widespread use as spatial light modu- 

lators, optical switchers and optical isolators. Recently it was realized that Bi:FG films with 
in-plane magnetization allow visualization and detection of magnetic fields • As basic 



sensors for magnetooptic imaging they have been used extensively in studies of magnetic 
flux in superconductors, domain formation in magnetic materials, currents in microelec- 



tronic circuits and recorded patterns in magnetic storage media |TJ-|18|. In spite of their 
successful applications, relatively few studies have been directed towards the understanding 
of the magnetooptical properties of Bi:FG films themselves. Of particular interest here is 
the behaviour and optimization of the Faraday rotation spectra. 

The dodecahedral sites in garnet crystals can accept several rare-earth ions, and each 
substitution seems to introduce its own special modification in the properties. All of the 
Bi^^, Ce^~^ and P6^+ have been shown to increase the Faraday rotation within certain wave- 
length intervals |ll8H2^. With Bi^^ a very strong Faraday rotation is obtained in the visible 
range, which together with the fact that the grown crystals are nearly transparent, make 
them by far the most frequently used material for magnetooptic imaging. For this appli- 
cation an important point is also that the films with in-plane magnetization are essentially 
free of magnetic domains. 

The understanding of the magnetooptical spectra in garnets is based on the knowledge 
of the spin orbit interaction in these materials. The spin orbit interaction generally causes a 
splitting of both the ground and excited states. However, in many cases it is possible to as- 
sume splitting of one of the states. Then only two situations are possible, (i) A paramagnetic 
transition, which is due to a orbital degenerate ground state with a singlet excited state, 
(ii) A diamagnetic transition, which is due to a orbital singlet ground state with a spin orbit 
split excited state. Since the paramagnetic and diamagnetic transitions have quite different 
lineshapes |]18[ , it should in principle be possible to distinguish between isolated diamagnetic 
and paramagnetic transitions, and therefore obtain a correct interpretation of the spectra. 
However, Bi:FG materials have several transitions within the UV and visible region, which 
complicates the interpretation considerably. This has lead to some controversy regarding 
the kind of transitions involved, and therefore how to reproduce the spectra theoretically 
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Simsa et al. measured the Faraday rotation of liquid phase epitaxy (LPE) grown Bi:FG 
films with small substitutions of bismuth . To interprete the spectra, they introduced one 
diamagnetic and several paramagnetic transitions. In this way they were able to correctly 
predict the behaviour of the spectra in the visible and near-infrared region. However, there 
were still some deviations between experimental and theoretical data when the bismuth 
substitution increased. More recently, Matsumoto et al. measured the Faraday rotation 
for Bi:FG films made by the gel method p^j27[] . Here reasonable agreement between the 
experimental data and the theoretical predictions was obtained with four paramagnetic 
transitions in addition to the four diamagnetic ones. However, the inclusion of paramagnetic 
transitions was not explained. Actually, Dionne and Allen questioned whether paramagnetic 
transitions can account for the contribution from Bi, since the strong superexchange field of 
the spin degeneracy exclude Zeeman splitting of the ground state, thus leaving a spin singlet 
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ground state [E4 



In this paper we report experiments and modelling of the Faraday rotation spectra in 
Bi:FG films suitable for magnetooptic imaging. To that end we have grown a series of 
Bi:FG films using the LPE technique, and characterised their chemical composition. The 
Faraday rotation spectra have been measured for photon energies between 1.77 and 2.53 eV, 
corresponding to wavelengths between 700 and 490 nm. It is shown that the Faraday rotation 
changes significantly with the amount of substituted gallium and bismuth. Furthermore, the 
comparison of experimental and theoretical data confirms that the magnetooptic response 
changes linearly with the bismuth substitution. 



II. EXPERIMENTAL 
A. Sample preparation 

Single crystal films of Bi:FG were grown by isothermal LPE from Bi20^/ PhO / B20^ flux 
onto (100) oriented gadolinium gallium garnet (GGG) substrates. The growth takes place 
while the 30 mm diameter substrate is dipped into the melt contained in a Pt crucible. 
During the growth, the parameters could be controlled to create low magnetic coercivity 
and in-plane magnetization in the garnet films. The thickness of the films was measured 
using a Scanning Electron Microscope (SEM), and their composition determined with an 
Electron MicroProbe (EMP). Thicknesses and compositions of four selected samples are 
listed in Table 1. For sample 1, 2 and 4 the film was grown on only one side of the substrate, 
whereas sample 3 has a 2 /im thick film on both sides of the substrate. 

From Table 1 one sees that the films can be represented by the following general formula; 
{Lu3_x-y-zYyBixPhz}[F e2-uaG0'Ua]{.F ^z-UiGau^)0i2, where {} indicates the dodecahedral 
site, [ ] the octahedral site, and () the tetrahedral site. Note that only the total gallium 
content u = Ua + Ud can be extracted from the EMP, and to determine Ua and Ud separately 
other techniques such as neutron spectroscopy must be applied. Usually, these techniques 
require bulk samples, and are not very useful for thin films. The films also contain small 
amounts of P6, typically of the order of z = 0.05. In the theoretical analysis of the Faraday 
spectra, we will neglect the contribution from Pb. 



B. Faraday spectra measurements 

The Faraday rotation as a function of wavelength at room temperature was measured 
using a setup consisting of a monochromatic light source, polarizer, two magnetic Helmholtz 
coils, sample, analyzer and detector. A schematic drawing is shown in Fig. |l]. The Helmholtz 
coils produce a magnetic field up to 1.1 kOe along the optical axis, which makes it possible 
to saturate the Faraday rotation of the sample. In this work all the samples were brought 
into saturation during the measurements. The Faraday rotation angle was measured as 0, 

F = , 

where dtot is the Faraday rotation measured for the combined film and substrate, dgub is the 
Faraday rotation of the bare substrate and t is the thickness of the film. The substrates used 
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in this study are 0.5 mm thick, and give only a small contribution of ~ 0.1 degrees. In the 
present study we found that interference effects due to finite thickness as well as the Faraday 
ellipticity could be neglected. The Faraday rotation has been measured across the useful 
visible wavelength range with an accuracy of 0.05 degrees. The wavelength was determined 
within ±1 nm. 

Shown in Figs. 2-4 is the observed Faraday rotation as a function of energy for samples 
1-4. Here is denoted by negative numbers since we follow the sign-convention of refs. 
P J2^J25[| . One sees that for sample 1-3, the magnitude of Bi;' increases monotonically with 
the photon energy. The largest rotation was found in sample 3 at 2.43 eV where Qp = 4.6 
deg/fim. Since sample 4 is relatively thin, we were able to measure the spectrum up to 
2.53 eV without reducing the accuracy of the measurement. This revealed a peak in the 
spectrum near 2.45 eV, where Bi? = 4.1 deg/fim. 



III. THEORY 

Our analysis of the Faraday rotation spectra is based on the molecular-orbital energy- 
level diagram of Dionne and Allen shown in Fig. |^. The basic mechanism behind 
the enhanced Faraday rotation is here the cooperative action of Fe^~^ ions with degenerate 
orbital terms that are split further by covalent interaction with Br'^ . In fact, substitution 
of bismuth results in an overlap between the 6p orbit of Bi^^, 2p orbit of and 3d orbit 
of Fe^~^. As discussed by Shinagawa and later by Dionne and Allen [^, this overlap 



results in an spin orbit splitting 2Aj associated with the tetrahedral and octahedral sites, 
see also Fig. ^. In the model presented here we ignore the contribution from transitions 
other than those introduced by bismuth substitution. In fact, LuIG and YIG have very 
similar Faraday rotation spectra with a positive peak located near 2.8 eV, but with very low 
Faraday rotation below 2.6 eV 0. This makes our assumption reasonable. 

In order to calculate the magnetooptic response, an expression which links the micro- 
scopic parameters to the Faraday rotation is needed. In magnetooptics, this is usually 
expressed through the complex permittivity tensor, which for Bi:FG has the form, 

/ eo ei \ 
V eo/ 

where ei = e'^ + je" represents the magnetooptical response and eo the non-magnetooptical 



part. In the visible region, the Faraday rotation can be written [24 



e.«£j.;. (2) 

where u is the the frequency, or equivalently the energy of the light, c is the speed of light in 
vacuum and n is the refractive index. Independent measurements show that it is reasonable 
to set n = 2.5 in the visible range PJ^. 

Assuming that the transitions are of electric dipole nature, the off-diagonal elements 
derived from first order time- dependent perturbation theory are given by |24 



4 



The inner sum is over the right- and left-hand circular polarized light, whereas the outer is 
over the optical transitions in the tetrahedral and octahedral sublattices. The prefactor (±) 
denotes a substraction. The a;j+ and tUj. represent the resonance energy for right and left- 
hand circular polarized light, respectively. Similarly, and /j_ are the respective oscillator 
strengths, while Fj is the half-linewidth of the transition. Furthermore, e and m are the 
electron charge and mass, respectively, whereas is the active ion density. 

To first order, it is reasonable to assume that N is directly proportional to the bismuth 



content x ||2^. This assumption has also been confirmed by experiments [§,0. Furthermore, 
it is known that the strong enhancement of Faraday rotation is caused by iron-pair transi- 
tions, involving both octahedral and tetrahedral transitions simultaneously |3y]. Therefore, 



iron dilution of either sublattice results in a reduction of the active ion density. For these 
reasons we propose that the active ion density can be written as, 

N = Noil - ua/3)il - uj2)x . (4) 

Nq is a constant, and may expected to be 1/3 of the density of rare-earth ions on the 
dodecahedral site, i.e. 1.3 x 10^^ cm'^/S. When x = 3, this interpretation implies that the 
dodecahedral site is fully occupied by bismuth. 

The oscillator strength for right- and left-hand circular polarized light are given by 



f^± = ^\<g^\X\e,±>? , (5) 

where h is Planck's constant, X is the electric dipole operator, and \gi >, \ei± > being the 
wavefunction of the ground state and excited state, respectively. Using that for the present 
case one has u!i± = tUj ± Aj, we obtain from eq. (|^), 

fi± = ^{cOi±A,) , (6) 



where /j is an effective oscillator strength |2^. This expression is based on neglecting the 
difference between |ei+ > and |ej_ >. 

In order to fully predict the Faraday spectra, one needs detailed knowledge about the 
quantum mechanical oscillator strengths. However, this is not easily accessible with so many 
ions involved, and we will instead treat them as adjustable parameters in the comparison 
with experimental data. With eq. (||) inserted in eqs. and (|^), one obtains 

nmc , UJi 

i=a,d ' 

juj^ + A^r-u'-T^ {u, - A.)^ - - F,^ ] 

1 [(a;, + A,)2 -uj^ + r^f + Auj'T^ [{u, - - + T^f + Au^^ J ' 

which represents our model for the Faraday rotation. 



IV. DISCUSSION 

To fit theoretical curves to the experimental data, the product N fi/x was chosen as free 
parameter. According to eq. (§) depends on the gallium and bismuth content. On the 
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other hand, fi is roughly independent on these two substitutions. Therefore, the product 
Nfi/x may be adjusted to any substitutions and transitions. The parameters Aj, oji and 
Fj were chosen as sample independent, and the values suggested by Dionne and Allen were 
used as a starting point in the fitting pSf. Table 2 presents the parameters we found to 



give the best fit between the theoretical curves and the experimental data. Note that our 
values for Ud-, and Ta differs slightly from those of Dionne and Allen. A reason for this is 
that the minimum Faraday rotation seems to be shifted slightly towards lower energies as 



compared to the data analysed in ref . . Furthermore our values for F^ and F^ are smaller 
to account for the narrow linewidth near 2.45 eV seen in Fig. ^. 

In Fig. 1^ and |^ the Faraday rotation calculated from eq. (0) is plotted as solid lines, 
and show an excellent agreement with experimental data. In Fig. |^ the dashed and dash- 
dotted lines represent the contribution to the total Faraday rotation from the tetrahedral 
and octahedral sites, respectively. Note that below 2.2 eV, the main contribution to Qp 
comes from the transition related to the octahedral site. 

Fig. ^ shows the theoretical curve (solid line) fitted to the experimental data for sample 
4. Again there is excellent agreement at small photon energies, but some deviation appears 
above 2.4 eV, where both the experimental and theoretical curves show a peaked behaviour. 
In fact the experimental data indicate a stronger peak than the theoretical prediction. There 
may be several reasons for such a discrepancy. First, eq. (^) represents a simplification for 
the present type of FGs. In particular, we have not considered the canting of the spins or 
the temperature dependence of A^. There is also a possibillity that Aj and Fj should be 
chosen differently to obtain a stronger peak around 2.45 eV. Finally, it must be noted that 
we have neglected the transitions associated with pure LuIG. 

We have in Fig. | also plotted theoretical predictions for the case x = 1.5 (dashed curve) 
and X = 2 (dash-dotted curve) assuming that gallium content is identical to that of sample 
4. When a; = 2 it seems possible to obtain a Faraday rotation exceeding 6 deg/fim at 2.3 
eV in this material. 

The distribution of gallium on the tetrahedral and octahedral sites have been examined 
by Czerlinsky and Scott et al. for various garnet-compositions with and without bismuth 
p9| , |30[] . Their results indicated that gallium will replace Fe^^ in both tetrahedral and 
octahedral sublattices, but with the strongest dilution of tetrahedrally coordinated Fe^^ . 
In fact, Czerlinsky found that for Y^Fe/^GaiOi2 the distribution is 0.9 Ga on tetrahedral 
sites and only 0.1 on octahedral sites. Although the samples investigated in this study have 
somewhat different composition, we will assume that the gallium distribution is the same. 
Then, based on eq. (^) the gallium dependence can be removed resulting in a Faraday 
rotation associated only with the Bi-content given by, 

^ {i-ud/m-uam ' ^ ^ 

where Ud = 0.9m and Ua = O.lu. In Fig. |^ the is displayed as a function of the bismuth 
substitution. The experimental data have the same notation as in the previous figures, and 
are shown for two different energies, 1.97 eV (630nm) and 2.3 eV (540nm), corresponding 
also to the theoretical curves drawn as solid and dashed lines, respectively. To obtain the 
theoretical curves we used N^fa = 2.77 x 10^^ cm~^ and Nofd = 6.30 x 10^^ cm~^, which 
are the average values of the four samples as calculated from Table 2 and the assumed 
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gallium distribution. The experimental points are located very close to the straight lines, 
thus indicating that our linear approximation is valid. However, the linear approximation 
may not hold for large x, since saturation phenomena could change the behaviour. In fact, 
experimental studies by Matsumoto et al. indicated breakdown of the linear approximation 
when x>2 



V. CONCLUSION 

We have synthesized single crystals of bismuth-substituted ferrite garnets by the liquid 
phase epitaxy method where GGG substrates are dipped into the flux. The growth param- 
eters are controlled to obtain films with in-plane magnetization and virtually no domain 
activity, which makes them ideally suited for magnetooptic imaging. The Faraday rotation 
spectra were measured across the visible range of wavelengths. To interprete the spectra 
we have presented a simple model based on the existence of two optical transitions of dia- 
magnetic character, one tetrahedral and one octahedral. We find excellent agreement with 
our experimental results within the visible range. It is shown that the Faraday rotation 
changes significantly with the amount of substituted gallium and bismuth. Furthermore, 
the experimental results confirm that the magnetooptic response changes linearly with the 
bismuth substitution. 

It would be of considerable interest to observe the behaviour at low and high tempera- 
tures, which will be a topic of future investigations. 
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FIGURES 



FIG. 1. Setup for measurement of the Faraday spectra. The analyzer is set at 45 degrees from 
extinction to maximize the signal. The polarize the light we use a sheet polarizer, whereas the 
analyzer is a Glan-Thompson prism. Together they have an extinction ratio of ~ 10~^. To measure 
the signal we used an ordinary silicon detector with 3x3 mm^ active area. 

FIG. 2. The Faraday rotation for sample 1. The squares are the experimental data whereas the 
solid line represents the theoretical fit. The dashed and dash-dotted lines represent the contribution 
from the tetrahedral and octahedral transitions, respectively. 

FIG. 3. The Faraday rotation for sample 2 (circles) and 3 (diamonds). The solid and dashed 
lines represent the theoretical fit for sample 2 and 3, respectively. 

FIG. 4. The Faraday rotation for sample 4. The dashed line represents a theoretical estimate 
of the Faraday rotation when the gallium substitution is identical to that of sample 4, but with a 
bismuth substitution of x = 1.5. The dash-dotted line is similar to the dashed line, but with x = 2, 
showing a peak value > 9 deg/ fim. 

FIG. 5. The basic molecular-orbital energy-level diagram. (-|-1) and (-1) represents right and 
left-hand circular polarization, respectively. Note that there are two transitions which influences 
the Faraday rotation; one associated with the tetrahedral site, the other with the octahedral site. 
The transitions are assumed to follow the selection rules for electric dipole transitions. 

FIG. 6. The Faraday rotation as a function of bismuth content. The experimental points have 
the same notation as in the previous figures, and are represented for two different energies, 1.97 
eV (630nm) and 2.3 eV (540nm). The lines are the theoretical curves for the same energies. 
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Sample nr. 


Lu 


Y 


Bi 


Pb 


Fe 


Ga 


t (//mj 


1 


1.619 


0.625 


0.683 


0.064 


3.842 


1.144 


7.5 


2 


2.141 





0.796 


0.065 


3.887 


1.088 


4.0 


3 


2.166 





0.874 


0.037 


3.937 


0.969 


4.0 


4 


2.186 





0.818 


0.062 


3.744 


1.170 


2.6 



Table 1: The thickness and chemical composition of the four samples. 
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site 






Nh/x 




A,; 








Sample 1 


Sample 2 


Sample 3 


Sample 4 








a 


1.72 * 10^3 


1.76* 10^3 


1.88 * 10^3 


1.69 * 10^3 


0.27 


3.15 


0.47 


d 


-3.90 * 10^2 


-4.00 * 10^2 


-4.26 * 10^2 


-3.85* 1022 


0.11 


2.51 


0.38 



Table 2: The parameters found to give the best fit between eq. (0) and the experimental 
data. N fi/x is measured in [cm"^], whereas Aj, uoi and Fj are given in [eV]. Note that the 
tetrahedral and octahedral sites are given different signs, since they contribute oppositely 
to the Faraday rotation. 
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